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Abstract 

The ratio of the depths of spectral lines is a powerful indicator of the effective 
temperature. The method based on this analysis is capable of discerning small 
temperature variations of individual stars. 

We apply this spectroscopic data analysis to three type of stars, namely an RS CVn 
type binary system, a young solar-type star and a Cepheid variable. We show 
that individual LDRs converted into temperature through calibration relations lead 
to rotational and pulsational modulation of the average surface temperature with 
amplitudes of 127 K, 48 K and 1466 K in the three types of stars, with average 
estimated errors of some tens Kelvin degrees. 

Keywords: stars: activity - stars: variability - stars: individual: HK Lac, HD 166, 77 Aql 
1. INTRODUCTION 

Data analysis in astronomy is crucial for the optimal extraction of scientific information. The rapid 
evolution of technologies related to telescopes and instrumentation needs to be complemented 
by more and more accurate techniques of astronomical data analysis. 

The determination of effective temperature plays a very important role in stellar physics. Indeed, 
this physical parameter is fundamental for locating stars in the HR diagram, for studies of gravity, 
mass, abundance and so on. Its measurement is not often as accurate as it is needed. For this 
reason it is important to find a suitable diagnostics of photospheric temperature and to develop 
an adequate data analysis software. 

The strength of absorption lines has been used for spectral classification since the first beginning 
of astrophysical spectroscopy, given its dependence on temperature and/or gravity, but it is not 
directly suitable for accurate temperature measurements. Although it is obvious that stellar spectra 
vary with temperature, the line strength depends on several parameters, making the interpretation 
more complicated than one might expect. The depth ratio of two absorption lines (LDR) with 
different excitation potential proves instead to be a very suitable diagnostics for this purpose, after 
an appropriate calibration has been done. It actually allows us to measure temperature variations 
with uncertainties of only a few Kelvin degrees ([8]). 

Our goal was to set up a method based onto LDRs to accurately monitor temperature variations 
that may occur during rotation of magnetically active stars or along the period of a pulsating star. 
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For this reason we used line pairs, in which one line is much more temperature-dependent than 
the other one or in which the depth variation with temperature has an opposite trend for the two 
lines. In fact, the ratio of the two line depth parameter is practically independent of rotational and 
macroturbolence broadening or spectral resolution since both lines are reshaped in the same way. 
Moreover, small differences in the spectrograph focus or scattered light do not affect LDRs, but, 
more importantly, differences in line strength coming from differences in chemical abundances 
largely disappear (for weak lines). 

We have performed a calibration of LDRs in absolute temperatures, taking into account the gravity 
effects, through a set of standard stars. Then, the individual LDRs of our targets have been 
converted into temperature through the aforementioned calibrations. 

The LDR-technique allowed us, at our medium resolution of R = 14 000, to measure the 
modulation of the average surface temperature in two classes of variable stars, i.e. spotted active 
stars and Cepheid pulsating variables, with sensitivity around 10 Kelvin degrees. 

By means of this data analysis, it is possible to derive the fractional spotted area (A rei ) and 
the spot temperature (T sp ) in the case of some spotted stars for which there are simultaneous 
spectroscopic and photometric observations. Moreover, for the study of pulsation properties, this 
spectroscopic analysis allows us to simultaneously record two fundamental parameters, namely 
temperature and radial velocity. 



2. OBSERVATIONS AND ANALYSIS 

We have acquired the spectra using the REOSC echelle spectrograph fed by the 91 -cm telescope 
at Catania Astrophysical Observatory - M. G. Fracastoro station (Serra La Nave, Mt. Etna). 
The spectral resolving power was about 14 000 and the thinned back-illuminated SITe CCD, 
1024x1024 pixels, allowed us to record five orders in each frame, spanning from about 5850 
to 6700 A. The average signal-to-noise ratio (S/N) at continuum in the spectral region of interest 
was 200-500 for the standard stars and 1 00-200 for the target stars. 

The data reduction was performed by using the echelle task of the I RAF 1 package following the 
usual steps: background subtraction, division by a flat field spectrum given by a halogen lamp, 
wavelength calibration using the emission lines of a thorium-argon lamp, aperture extraction and 
continuum fitting with a low order polynomial. Additional information about data reductions are 
given in [2]. 

3. MEASUREMENTS OF LINE-DEPTH RATIO 

The studies of temperature variations on active stars made by Gray and his collaborators 
are principally based on the lines pair Vi A6251.83 A-Fei A6252.57 A. Within the spectral 
range covered by our echelle frames there are several pairs of line suitable for temperature 
determination. The line pairs for each ratio have been chosen to be close one another in order to 
minimize errors in setting the continuum and in a region of low contamination from telluric lines. 

In Table 1 we list the wavelengths (Ai, A 2 ), the lower excitation potentials (in eV) of the selected 
lines ([13]) and the working temperature range of each LDR-r eff calibration (see Sect. 4 for 
description of the calibration). We have considered metal lines because they are the most 
sensitive to the temperature. Altogether we have identified 30 spectral lines forming 20 pairs 
suitable for line-depth ratios. These lines were identified through the solar spectrum atlas ([13]) 
choosing unblended lines. The only exceptions are four lines, each composed of two very close 
lines with the same temperature dependence, but not resolved at our resolution and behaving as 
single lines. 

1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of the 
Universities for Research in Astronomy, inc. (AURA) under cooperative agreement with the National Science Foundation. 
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TABLE 1: Spectral lines for developing calibrations using line-depth ratios r = d\ 1 /d\ 2 , with c/ line depth. 
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At our resolution, the center of a line profile and therefore the true depth may not concide with an 
observed data point. So we have developed an IDL {Interactive Data Language) routine to define 
the line depth. The lowest five points in the core of each measured line were fitted with a cubic 
spline and the minimum of this cubic polynomial was taken as the line depth (see Fig. 1). Writing 
the line depth das 

d= ^ = 1 - TT , (1) 

where S c and S b are the signals in ADU (Analog to Digital Units) or in photons of continuum and 
bottom of the line, respectively; the fractional error on dean be expressed as 

g d _^ _i-d rr~± (2 ) 

d d d V S b S c ' 

Assuming that signals in the continuum and in the bottom of the line follow a poissonian statistics, 
as in the case of a linear detector like a CCD, and considering, as an example, a line depth d = 
0.40 and values of 45000 ADU and 28000 ADU for the continuum and the bottom of the line of 
a generic raw spectrum (see Fig. 2), given the CCD conversion factor 2.5 e~/ADU we find this 
typical relative error: 



-i « 0.72% . 
d 

As we can see below, this can ensure temperature accuracy of 10-20 K. 
4. SPECTROSCOPIC DATA ANALYSIS 



The calibration of LDRs to temperature has been made by means of spectra of 43 standard stars 
with different spectral type (F6-M4) and luminosity class (V— 1 1) acquired from 2000 to 2002. Since 
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FIGURE 1 : Example of normalized spectrum of the standard star a Ari in the 6250 A region. The red 
line is the polynomial fit on the lowest five points in the core of the Fei A6252.57 A line. The green 
segment on this line has the signals in the continuum and in the bottom of the line as extremes. 
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FIGURE 3: An example of effective temperature as a function of LDR. The solid 
line represents a fit to evolved stars LDR, while the dashed line is the polynomial 
fit to MS LDR. 

effective temperatures are available for few of our calibration stars, we have adopted effective 
temperature derived from the empirical relation (B-V)-T eff proposed by [7]. Then, the calibration 
curves for main sequence (MS) and giant stars have been made separately. 

In Fig. 3 we show the A6275 VI/A6270 Fel LDR- 7"^ calibration as an example. The two different 
behaviours, for MS and evolved stars, are evident. The gravity dependence of LDRs has been 
deduced from the residuals of the standard star LDR compared to the fits of Fig. 3. These 
differences correlate very well with the gravity indicator A/W v (Fig. 4), which is the absolute 
magnitude difference in the HR diagram compared to the ZAMS. The best-fit in Fig. 4 has been 
used to correct the measured LDRs for gravity. A gravity-corrected LDR-r eff calibration for A6275 
VI/A6270 Fel ratio is shown in Fig. 5 for MS stars. Data points for stars of different luminosity 
classes mix well together, allowing a unique correlation, alternatively applicable to MS and giant 
stars. The temperature sensitivity of each line-depth ratio derived from the slope of the polynomial 
fit is of the order of 10-20 K for 1% LDR variation. 

5. LINE-DEPTH RATIO METHOD APPLIED TO THREE TYPES OF STARS 

LDR or temperature variations can occur when hotter or cooler surface features are carried across 
the disk of an active star as the star rotates or during natural oscillations of a Cepheid star. 

On the basis of these considerations we have made observations on three type of stars in order 
to investigate the applicability of the LDR-analysis. 

5.1. RESULTS ON A RS CVn BINARY 

As an example of the application of this analysis technique, Fig. 6 shows the temperature 
modulation of the slow-rotating single-lined active star HK Lac (HD 209813, KOMI) during its 
rotational period of P rot = 24 d .4284 ([6]). The spectra have been acquired from October 2000 to 
January 2001 and the initial epoch is taken from [15]. The figure displays an average of six LDR- 
curves converted into temperature curves adopting the calibration for giant stars. The amplitude 
of the < T eff > curve we measured is about 130 K and the corresponding errors are tipically of a 
few Kelvin degrees. This temperature variation is the result of the different visibility of cool spots 
along the rotational period. The upper box of Fig. 6 also showns the synthetic curve, obtained with 
a two-spots model (lower panel), superimposed on the data points. 
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FIGURE 4: Residuals of LDRs compared to the polynomial fit to MS data 
of Fig. 3 versus Al\/kj. The intercept and the slope of the regression line 
are also reported. An analogous regression line has been obtained for 
LDRs compared to giant's polynomial fit. 
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FIGURE 5: Effective temperature as a function of corrected LDR for MS stars. 
The solid line represents a fit to all data, while the blue dashed line is the 
polynomial fit to LDR data corrected for giant-star calibration. A similar fit of 
T e fF has been obtained as a function of corrected LDR for giant stars. 
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FIGURE 6: Observed (dots) and synthetic (full line) average effective temperature variation of HK Lac for the 
two-spots distribution shown in the lower panel at four specific phases. The minimum of temperature curve 
is around the phase 4> = P .00 when the two spots are both visible on the disk of the star. 
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FIGURE 7: Upper panel: temperature curves of HD 166 obtained from each 
LDRs. Lower panel: effective temperature as a function of the rotational phase. 



5.2. RESULTS ON A SOLAR-LIKE STAR 

We are also studying the rotational modulation of the temperature through LDR in young solar- 
type stars. They exhibit intermediate activity levels between that typical of RS CVn binaries and 
of the Sun. 

Fig. 7 reports preliminary results about the temperature variation of the slow-rotating single main 
sequence star HD 166 (K0V) during its rotational period of P rot = 6 d .23 2 . The observations have 
been made in 2000 from October 16th to November 16th. HD 166 shows a clear rotational 
modulation of the average temperature, obtained from 8 LDRs, with an amplitude of about 50 K. 

5.3. RESULTS ON A CEPHEID STAR 

We have also applied the LDR-method to the Cepheid star r) Aql (F6lb-G4lb, P pu | S = 7 d .1 76641 3 ). 
The observations have been performed in 2002 between July 22nd and October 6th. 

In Fig. 8 we plot the light curves obtained by [1] together with those of [9]. The difference between 
the two data sets does not exceed 0.01-0.02 mag, which is within the explicitation in matching 
different data sets. 

In the same figure we show the effective temperature obtained by us (dots) and by [1 0] (diamonds) 
by means of the LDR-method and that derived by [4] (asterisks) using the independence from the 
excitation potential of Fei lines of iron abundance calculated using Kurucz LTE, plane-parallel, 
atomic and molecular line-blanketing, flux-conserving stellar models. The agreement between 
these three data sets is good, the difference of values being at the same phase of only few per 
cent of the total amplitude (A < T eff >~ 1500 K). 

2 We have used the mean period and the initial julian day given by [5]. 
3 We have used the period and the epoch of brightest Vgiven by [12]. 
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FIGURE 8: Upper panel: light curves obtained by [1] and [9]. Middle 
panel: our, [4]'s and [10]'s average effective temperature as a function of 
pulsational phase. Lower panel: radial velocity curves from our and [9]'s 
observations. 



From the same spectra we derived the simultaneous radial velocity curve, displayed in Fig. 8 
together with that of [9]. Also in this case, the two sets of points are almost superimposed and 
the maximum difference compared to our curve amplitude of about 40 km/s is 5%, i.e. within the 
typical accuracy of RV measurements. 

We can also note that the radial velocity curve of r\ Aql is almost a mirror image of the light 
curve, so that the luminosity maximum corresponds to the highest velocity of recession. These 
two curves are lightly shifted, as predicted by the pulsation theory ([11]). 

6. CONCLUSIONS 

The line-depth ratio method proved to be very effective for measuring the surface temperature in 
several classes of stars. We have shown that echelle spectra, covering a large spectral band and 
allowing the use of many line pairs, enable us to enhance the precision of effective temperature 
measurements to few degrees. 

Well-defined modulations of the average temperature have been detected in HK Lac, HD 166 and 
j] Aql. 
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HK Lac displays a clear rotational modulation of its surface temperature as derived from several 
LDRs (see [2] for major details). From the combined analysis of contemporaneous temperature 
and light curves, by applying a spot-model, we can obtain sole solutions for the spot temperature 
T sp and the relative area coverage A re \ with an accuracy of few per cent ([3]). Spot temperatures 
we derive are closer to those typical of sunspot penumbra rather than umbra. 

For HD 166, our analysis represents a snapshot of its active regions and more observations are 
required to study any time dependence. It is important to make observations of young solar-type 
stars because their study provides insight into the behavior of the solar activity and its effect on 
the planets during the first few hundred million years of solar system history. 

We have also derived a strong modulation of 77 Aql temperature during its pulsational phase. 
Determining accurate values of T eff in Cepheid stars along their pulsational cycle is of paramount 
importance. For example, precise temperatures provide a reliable localization of variable stars in 
the Cepheid instability strip, whose width is only about 650 K. Moreover, the combined analysis 
of temperature, radial velocity and light variations coming from spectroscopic and photometric 
data allows us to derive accurate values of Cepheid absolute magnitudes. This parameter is 
of fundamental assistance for distance determination, since the Cepheids are considered the 
primary standard candles together with the RR Lyrae variables. 
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